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Abstract 

Various methods of constructing an orthonomal set out of a given set of lin- 
early independent vectors are discussed. Particular attention is paid to the 
Gram-Schmidt and the Schweinler-Wigner orthogonalization procedures. A 
new orthogonalization procedure which, like the Schweinler- Wigner proce- 
dure, is democratic and is endowed with an extremal property is suggested. 
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Constructing an orthonormal set out of a given set of linearly independent Vectors is 
an age old problem. Among the many possible orthogonalization procedures, the two algo- 
rithmic procedures that have been extensively discussed and used in the literature are (a) 
the familiar Gram-Schmidt procedure [1] and (b) a procedure which is referred to as the 
Schweinler -Wigner procedure [2], particularly in the wavelet literature [3]. ( This method is 
known among the chemists as the Lowdin orthogonalization procedure [4]. Mathematicians 
attribute it to Poincare. Schweinler and Wigner themselves trace its origin to a work of 
Landshoff [5]. Eschewing the question of historically correct attribution, we shall continue 
to refer to it as the Schweinler- Wigner procedure ) An intrinsic difference between the two 
procedures is that while the Gram-Schmidt procedure, by its very nature, requires one to 
select the linearly independent vectors sequentially, the Schweinler- Wigner procedures treats 
all the members of the set of linearly independent vectors democratically. The significance of 
the work of Schweinler and Wigner lies not in introducing a new orthogonalization method 
- the method was already known but rather in introducing a positive quantity m, to be 
defined shortly, which discriminates between various orthogonalization procedures. They 
showed that m is a maximum for the Schweinler- Wigner basis. In this letter we pose and 
answer the question as to what is the orthogonalization procedure which minimizes m. This 
new orthogonalization procedure, like the Schweinler- Wigner procedure, also turns out to 
be completely democratic in that it treats all the linearly independent vectors on the same 
footing. The quantity m was introduced by Schweinler and Wigner in a some what ad-hoc 
manner. We reformulate their procedure in a way so as that the quantity m appears in a 
natural way and can be useful in a wider context than that for which it was introduced. 
In particular, this reformulation enables us to quantify the notion of an orthonormal basis 
which brings any Hermitian operator in to a maximally off-diagonal form. 

Let • • • , v N denote a set of N linearly independent vectors. Let M denote the associ- 
ated Gram matrix :Mij = (vi,Vj). M is a positive definite Hermitian matrix. Define 

z = vS , (1) 
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where S is an invertible matrix. Then 

(zh z j) — (S^MS)ij . (2) 

Requiring that z be an orthonormal basis amounts to requiring that 

S^MS = I i.e. M- 1 = SS^ . (3) 

Each such S defines an orthogonalization procedure. Two standard choices of S are 
[1] Schweinler-Wigner Procedure 

This procedure corresponds to the choice 

S = UP- 1 ' 2 tf , (4) 

where U is the matrix which brings M to a diagonal form P 

U^MU = P . (5) 

With this choice of S, which corresponds to taking the Hermitian square root of the matrix 
M, one has 

z = vUP- 1 ' 2 tf . (6) 

On defining w = zU, one obtains the Schweinler-Wigner basis 

w = vUP- 1/2 . (7) 

[2] Gram-Schmidt Orthogonalization Procedure 

In this procedure S is chosen to be an upper triangular matrix T satisfying 

and the Gram-Schmidt basis is given by 

y = vT (9) 



The two orthonormal bases w and y discussed above are related to each other by the 
following unitary transformation 

y = wV (1) , (10) 

where is given by 

y(l) = pV*rf T . (11) 
Schweinler and Wigner introduced a quantity m(z) as follows 

(*) = E fe I (^) l 2 V , (12) 

k \ I ) 



where z is any orthonormal basis. They further showed that m(z) attains its maximum 
value Tr(M 2 ) for z = w. 

m max =Tr{M 2 )=m{w) . (13) 
For any other basis z, related to w by a unitary transformation V 

z = wV , (14) 

the value of m(z) is given by 

™ ( z) = E((^)J • (15) 



( ( v j pv) ' 2 

k 

In particular, for the Gram-Schmidt basis y, one finds that 



"(y) = E(WJ' ■ (16) 

k 



A natural question to ask is as to what is the orthonormal basis which minimizes m(z). 
On applying Cauchy inequality to (15), it follows that 



m(z) = £ ((V'PV)J 

k 



>l(TrM) 2 . (17) 



Equality sign holds if and only if (V^PV)kk = c independent of k i.e 

1 



ki 



N 



(Pi + • • • + P N ) 



(18) 



This requires that 



V k i | = — 



(19) 



for all k and /. The matrix elements of the unitary matrix satisfying the above equation 
are thus given by 



(2) 



kl 



cxp 



2m(k - - 1] 



N 



(20) 



The matrix is thus just the character table of the cyclic group Cn- Thus the basis x, for 
which m(z) attains its minimum value m m i n = (l/N)(TrM) 2 , is related to the Schweinler- 
Wigner basis w as follows 



Xl 



E ex P 



2m(k-l)(l-l) 
N 



(21) 



This basis also treats all the linearly independent vectors v democratically like the 
Schweinler-Wigner basis. 

The quantity m(z) appears to have been introduced by Schweinler and Wigner in a rather 
ad-hoc way. At least no particular motivation for introducing it appears in their work. We 
now reformulate their work in a way that this quantity appears naturally. Consider the 
Hermitian operator 



In an arbitrary orthonormal basis z one can write 



(22) 



Tr(M 2 ) = J2 I (zi,Mz m ) | 2 + EI (zuMzi) 



(23) 



Im 



The second term on the rhs is easily seen to be the same as m(z) in (12). From this 
perspective it is immediately obvious that the basis which maximizes m(z) is the one in 



which M is diagonal. This is just the Schweinler-Wigner basis as can also be directly 
verified. Thus the Schweinler-Wigner basis is simply the eigenbasis of the operator M. and 
if the eigenvalues of M, are all distinct then this basis is essentially unique. Further, since 
Tr(M 2 ) is independent of the choice of basis, it is clear from (23) that the basis which 
minimizes m(z) maximizes 

n(z) = ^2\^i,Mz m )\ 2 . (24) 

Im 

The quantity n(z) therefore provides a quantitative measure of the off-diagonality of the 
operator M. in the z basis. The new orthonormal basis proposed in this work is thus the 
one in which Ai is maximally off-diagonal. 
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